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1. Introduction and Summary 

Cosmology acts as a test bed for theories beyond the Standard Model, including theories 
for quantum gravity and string theory. The current observational data strongly favors 
primordial inflation, close to the grand unified scale Q. Inflation is a dynamical solution 
to the flatness, horizon and homogeneity problems Q. In a generic model, such as one 
driven by a scalar field, it also stretches metric fluctuations outside the Hubble radius with 
a scale invariant power spectrum [^. 

In spite of the great successes of inflation, it does not address one of the most important 
aspects of the Big Bang cosmology: For any equation of state obeying the strong energy 
condition p > —p/3, regardless of the geometry (flat, open, closed) of the universe, the scale 
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factor of the universe in a Friedmann Robertson Walker (FRW) metric vanishes at t = 0, 
and the matter density diverges. In fact all the curvature invariants, such as R, OR,..., 
become singular. This is the reason why it is called the Big Bang singularity problem. 

Although inflation requires p < — p/3, it does not alleviate the Big Bang singularity 
problem, rather it pushes the singularity backwards in time. Many authors have pondered 
on this debatable issue, whether inflation is past eternal or not ||], and the conclusion is 
that it is not, at least in the context of Einstein gravity as long as the average expansion 
rate in the past is greater than zero, i.e. Hav > p. The fluctuations grow as the universe 
approaches the singularity, and the standard singularity theorems due to Hawking and 
Penrose hold, which inevitably leads to a collapse in FRW geometry as long as the energy 
density is positive Q (see also P). 

There have been other attempts to circumvent the Big Bang singularity problem due 
to anisotropic stresses, self regenerating universes (during inflation), quantum cosmology, 
etc. (see [^) but none has successfully resolved the issue of space-like singularity, espe- 
cially in the context of a flat universe^. Even string theory has yet to address this issue 
comprehensively: Several toy model constructions have been attempted, but mostly they 
encounter some pathologies, such as closed time-like curves, quantum instabilities, presence 
of ghosts, negative-tension branes, singular bounce, etc. (see for instance @)^. 

In this paper we will rather seek a phenomenological solution, but without any patholo- 
gies. We will seek a bouncing solution within 3-1-1 dimensions where the scale factor of a 
flat, homogeneous and isotropic metric 

ds'^ = -dt^ + a'^it) [dx^ + dy^ + dz^] (1.1) 

undergoes a non-singular bounce, a course of initial contraction to a minimal radius, then 
a subsequent phase of expansion to enter into the hot Big Bang era. 

We will advocate higher-derivative corrections to the Einstein-Hilbert action that are 
both ghost free and asymptotically free to find sensible non-singular bouncing solutions 
in the presence of a fluid with non-vanishing pressure. Earlier attempts to find such an 



action were not very successful |11] (see also p^ , |ll3[| ). Below we briefly explain the kind of 



actions that we consider and why we believe they can circumvent the problems encountered 



in |11], [12|, and others: 



Stringy motivation: 

The main theoretical motivation comes from string theory, which suggests higher- 
derivative corrections to the Einstein-Hilbert action. There such corrections appear 
already classically (i.e., at the tree level), but we do not preclude theories where such 



^In the context of a closed universe where the curvature term acts as a "source" for negative energy 
density in the Hubble equation, one can obtain bouncing solutions but this additionally requires a phase 
of inflation, to address the "flatness problem". Although this is a perfectly viable scenario, the motivations 
for a bounce is a little bit lost, as it is also supposed to be an alternative to inflation. 

^In jlo| the authors considered several brane world models. One such example is where a probe anti- 
brane circles around a stack of source branes giving rise to an induced 3+1 dimensional metric. The bulk 
metric uausally has a horizon and also singularities, but the induced metric can undergo a non-singular 
bounce. One of the concerns of these models is the issue of stability and lack of full control of the system. 
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corrections (or strings themselves) appear at the loop level or even non-perturbatively 
(e.g., as expected in the expansion in some Yang- Mills theory From string 



field theory [15| (either light-cone or covariant) the form of the higher-derivative 
modification can be seen to be Gaussian: There are factors appearing in all 
vertices (e.g., (e'-'i;^*)^), which can be transplanted to kinetic terms by field redefinitions 
((/) — > e~^(f)). The nonperturbative gravity actions that we consider will be inspired 
by such kinetic terms suggested by string theory. 

• Ghost free action: 

Previously, non-singular bouncing solutions have been obtained by introducing states 
carrying negative energy compensating the positive energy contribution, by introduc- 
ing new fields or higher derivatives (see for instance [0]). In general, such theories 
are plagued by ghosts - states that violate unitarity. Alternatively, ghosts carry neg- 
ative energy which inevitably leads to vacuum instability |T7[| . Thus, one condition 
we impose is the absence of ghosts. Note that the Gaussian factors of string theory 
do not introduce new states (ghosts or otherwise) because they have neither poles nor 
zeroes at finite momentum (which is also why they are allowed in field redefinitions) . 

• Asymptotic freedom: 

This is required (in addition to renormalizability) in ordinary field theory for phe- 
nomenological reasons (scaling in processes such as deep inelastic scattering), as well 
as to avoid certain problems in the nonperturbative definition of the theory. (This 
can be seen, e.g., in lattice quantum chromodynamics when defining a continuum 
limit. In the case of an effective theory for gravity, this condition takes the form 
of "asymptotic safety" [|^, a more general statement of a well-defined ultraviolet 
limit.) In gravity this feature provides a solution to the singularity problem; here 
we address the Big Bang singularity. (Black hole singularities were considered in 
[p^].) Asymptotic freedom implies that gravity becomes weak at short distances or 
high energies, enabling pressure to counteract the gravitational attraction, thereby 
preventing a collapse. Again the Gaussian factors are more than enough to produce 
such behavior. (The fact that this fall-off is so much stronger than normal asymptotic 
freedom is why string theory cannot describe parton behavior as a model of hadrons 
for QCD.) 

We will find that if one wants to have both a ghost and asymptotically free theory of 
gravity, one has little choice but to look into gravity actions that are non-polynomial in 
derivatives. To see this let us start by considering a simple generalization of the Einstein- 
Hilbert action, the so called fourth-order gravity: 



S = j (fx V^F{R) , (1.2) 

with 

F{R) = R + cqR^ + boC^ . (1.3) 
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where C refers to the Weyl tensor. This action, in general, is asymptotically free. In fact, 
the improved ultraviolet behavior even makes it renormalizable [^]. Unfortunately this 
theory also has a (Weyl) ghost (see also [|l3|), unless ho = 0, in which case the nice 
property of asymptotic freedom is also lost. Similar conclusions can also be verified from 
the analysis done by |^0[ for more general actions involving the Ricci and the Riemann 
tensors, when F = F{R,R^^,R^''' ^R^^p^m'^P"). 

In the most general case one can consider actions involving the scalar curvature, the 
Ricci and Weyl tensors, as well as with arbitrary higher derivatives on these curvatures. If 
the number of derivatives (including those in the curvatures) is finite, the kinetic operator 
will be polynomial, and by the usual arguments will (1) contain a single pole, (2) contain 
a ghost, or (3) not fall off at large momenta faster than a single pole (i.e., for the case of 
gravity, not be asymptotically free). This argument can be applied separately to the spin 
and spin 2 pieces of the metric; only the spin piece contributes to cosmology, but the 
spin 2 piece contributes to black holes. In section |^ we discuss this for the spin piece 
(scale factor) , which means the scalar curvature is the only non- vanishing one (conformally 
flat space). 

Therefore, in this paper we mainly focus on actions which contain arbitrary powers of 
higher derivative terms: The simplest such action is of the form 

oo 

F(i?) = ii + ^c„i?"+2. (1.4) 

n=0 

Unfortunately, even such actions cannot simultaneously be ghost and asymptotically free, 
which is easily seen from its correspondence with scalar-tensor/Brans-Dicke kind of theo- 
ries (see for instance [^]). Nevertheless, such theories are cosmologically interesting alter- 
natives to Einstein-Hilbert gravity, and can be tested by the solar system constraints or 
through spinning objects such as pulsars [p^], or even cosmologically p^]. Such theories (or 



their scalar-tensor analogues) have been studied in the context of inflation 24], creation 
of the universe through an instanton ||25[, understanding reheating after inflation [ |26| , and 
more recently, understanding the origin of dark energy (see for instance [^] ) . In section ^ 
and appendix we review their dynamics in the context of realizing a past asymptotically 
de Sitter universe and also point out why it doesn't constitute a resolution of the Big Bang 
singularity, although it may be relevant for inflation. 

The next simplest non-perturbative action is of the form^: 

oo 

F{R) = R + ^CnRW'R (1.5) 

n=0 

We will show that such actions can indeed give rise to a ghost and asymptotically free 
theory of gravity, and this is one of the key results in the paper. Moreover, we will also 
obtain exact bouncing cosmological solutions for such actions, thereby addressing the Big 
Bang singularity problem, see section 0. We note in passing that (|1.5|) should be treated as 



Such actions can also be motivated by symmetry considerations, such as local (Weyl) scale transforma- 



tion, and in section 3.3 we provide such an example. 
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an "effective action" , and thus the fact that it contains an infinite number of derivatives, 
and is therefore non-local, is only to be expected. 

We also study the Newtonian potential for the action given in ( |1.5D ; we show that at 
large distances and at late times contributions from the higher derivative terms become 
negligible and the theory is effectively given by the Einstein-Hilbert action. Therefore, 
the theory possesses Minkowski space as a low energy vacuum solution, recovering the 
Newtonian potential with 1/r dependence, where r is the length separation. On the other 
hand at small distances the higher derivative terms play, more and more, an important role, 
thereby indicating the importance of asymptotic freedom near the bounce: see section |5|. 

The paper is organized as follows: In section |^ we mostly review difficulties in obtaining 
a singularity free universe. Partially for simplicity, we restrict ourselves to a spatially flat 
universe. Specifically in the context of F{R) theories we find that in the absence of ghosts, 
no satisfactory singularity free evolution is possible, corroborating our earlier arguments 
based on the lack of asymptotic freedom. We do note however that a small class of actions 
that grow at most as i?^ as — > oo can provide a past asymptotically de Sitter solution 
similar to what was found previously in for (1.3) and 60 = 0. In the following sections 
therefore we focus on actions of the form ([1.5|). 

In section ^ we explain when actions of the form ( |1.5| ) can be both ghost and asymp- 
totically free. We also obtain the low energy Newtonian limit. In section ^, we obtain exact 
and approximate bouncing solutions, and in the next section ^ we show how a transition 
to normal FRW evolution occurs before and after the bounce. 



2. F{R) Cosmology 



In the field theoretic framework a resolution to the Big Bang singularity would imply 
a singularity free history from t ^ —00 till today. Here t denotes the proper time. So the 
picture is that the universe simply exists and one can go as far back in time as one wants 
without encountering any singularity. A priori, three different evolutions are possible: 

I. The universe expanded monotonically. 

II. The universe first contracted and then expanded (in other words, a "bounce" solution). 

III. The universe undergoes periodic phases of contraction and expansion, or a cyclic 
universe scenario. 



Since III requires multiple bounces and turnarounds, it also obviously entails II, and 
therefore we will in this paper only focus on determining whether I or II are possible. 

In this section we study cosmological evolutions with F{R) theories. We find that 
only a very restrictive class of such functions can provide a singularity free history of our 
universe. To show this we proceed in the following steps: First we explain the connection 
between F{R) theories given by ( [1.4| ) and scalar-tensor theories |28|. Next, we study 
scalar-tensor theories and find out for what kind of potentials one can have singularity free 
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cosmologies. Finally, we translate back the scalar-tensor action to find which F{R) actions 
can provide a singularity-free evolution this way, but we also point out why they cannot 
address the Big Bang singularity problem comprehensively. 

2.1 Conformal Transformation and Singularities 



We start by noting that an action of the form (^]^) is equivalent to a Brans-Dicke 
theory [^] of the form 

Sbd = J d^x ^[^R - Vm , (2.1) 
where V{^) is defined through a Legendre transform 

V{^) = ^R{^)- F{R{^)), (2.2) 
with _ 

(.3) 

defining R = R{^). Note that for a ghost free action $ is always positive; it was shown in 
pO| ], that when F'{R) is negative, one is plagued with ghosts. This fact then allows us to 
perform a conformal transformation 

= ^9t.u = eV^5<^<7M. , (2.4) 
so that the action is that of a scalar-tensor theory 



where 



R! - (d'cPY - 2 



Ml , 



(2.5) 



y(0) = ie~V §<^i?(0) - ie" V "^'^F{R{(t>)) . (2.6) 
Thus, at least classically, ghost-free F{R) theories are physically equivalent to the scalar- 



tensor theories given by (2.5). 

So, let us now ask when can one have a singularity free evolution in the context of 
scalar-tensor theories. It is known already that such theories with ordinary (non-ghost like) 
kinetic terms as in ( p.5D do not permit a bouncing cosmology, as the sum total contribution 
to the energy densities from the scalar and the gravity sector remains positive. Therefore 
one does not obtain a bouncing universe^, i.e., situations summarized in II, III. 

The only option left is I, i.e., to find a monotonically expanding solution. Indeed, 
such solutions do exist within scalar-tensor theories (see for instance ^]). However, 
in appendix ^ we provide a general argument to show that invariably such a scenario 
corresponds to an asymptotically de Sitter universe in the past which is known to be 



*There is an intriguing possibility that although our universe may now be in a ghost free state, F'{R) > 0, 
it could have passed through an earlier phase when F'{R) < 0. Can such a "transient ghost" degree of 
freedom cause a bounce? 
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geodesically incomplete @ and therefore fails to address the Big Bang singularity in a 
comprehensive manner. It may still be interesting from the point of view of inflation. In 
appendix^, we show that: only a subclass of F{R) theory can provide such an inflationary 
solution, only F{R) 's which grow at most as R'^ as R ^ oo. 

3. Theories with F{R, OR, . . .) 



Since ghost-free F{R) type actions do not satisfactorily resolve the cosmological Big 
Bang singularity, we therefore look at more general actions of type ( |1.5D including OR type 
of terms in the Lagrangian. 

3.1 Ghosts 

Although it is commonly believed that higher derivative theories contain ghosts, this is 
not always true. For example, the F(R) actions considered in the previous section do not 
contain any ghosts. In Lorentz covariant gauges there appears to be a scalar ghost, along 
with the usual ghostlike (negative-metric) time components of the metric. However, this 
ghost already appears in the usual Einstein-Hilbert action: It is the scale factor of the 
metric. It disappears in the lightcone gauge, along with the longitudinal components of 
the metric. The addition of the R^ term then contributes a scalar mode of opposite sign, 
which is thus physical: the right sign, and not a gauge artifact. The same procedure would 
not be possible with additional higher derivative contributions to the propagator, which 
would give spin or 2 with the wrong sign in any gauge. A simple way to analyze this 
effect for spin modes is to study the equivalent scalar-tensor theory where all the higher 
derivative terms reside in the scalar sector, and one has to worry about ghosts in the scalar 
sector only^. 



We start with the quadratic action given by (1^). This action is equivalent to a higher 
derivative scalar-tensor action, given by 



S 



d^Xy 
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^R + 7PJ2 Ci^V - {^(^ - 1) - co^^} 



(3.1) 



The easiest way to see that ( p.l| ) corresponds to (|l.5|) is by looking at the field equation 
for $: 

6 S 

— = 0^V = ^- (3.2) 

Then by substituting ip in (|3.1| ) one recovers ( |1.5D . We now perform a conformal transfor- 
mation = $^/^e^™, and note that dip = O'-ip + 0{(f>'^ , 4>^p,'ip'^) where we have defined 
^ = e't'r^l-\-(j)-\- 0{4>^). Then up to quadratic terms we find 



S w d^x 



R' + I^DV + ^ qD'V - - co^^} 



(3.3) 



''Alternatively, one can restrict oneself to conformal metrics, and just compute the propagator for the 
conformal scale factor. 
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To look for ghosts we have to find the propagators. The field equations for (p and ij: are 
given by: 

ip = 'iU(t), (3.4) 

" oo 

^qDV + co^ . (3.5) 
. 1 

Substituting il) from (3^) in ( p. 5] ) we find: 

(i-6^c^n'+M (/> = r(n)(/) = o. (3.6) 



From (3^) we can easily read off' the scalar propagator to be: 

1 

r(-p2) 



Gip") - T^r-JTv (3-7) 



Let us consider the contrasting cases. 

• Examples of ghosts: 

First we assume that F is a finite power series. In this case one can always write F 
as: 

F( V) ~ ip" + mDip" + ml)... (p2 + ml) . (3.8) 

In order for the theory to be non-tachyonic, all the mf have to be positive and real. 
Moreover, if there are at least two discrete single poles (say mi ^ W2), then at least 
one of them is ghost like (one of the residues has to be negative). A double pole can 
be represented as the convergence of two simple poles with opposite residues. Similar 
arguments follow for higher order poles. 

The above argument can be generalized to any continuous F(— p^). Suppose F(— p^) 
has at least two distinct zeroes. Say = —m\ and = —m^ are two adjacent ones, 
with < m^. Then it is of the form 

F( V) = (/ + mi){p^ + m2)/( V) • (3.9) 

Since the zeroes considered above are adjacent there are no more zeroes in /(— p^) 
between — m| < p^ < — m^. Thus the sign of f{—p'^) has not changed in this range. 
It then follows that the residue at p'^ = —m\ and p^ = have different signs. 

Therefore again there is a ghost. 

• Requirement for a ghost free case: 

For a polynomial inverse propagator (^^) the only case when the propagator is ghost 
free is 

F(n) = n-m^, (3.10) 

which corresponds to the R? type of action that we mentioned in the previous section 
(also, see appendices 0,^. More generally, one can at most have a single zero in 
F(V). 



-8- 



As we argued before, stringy redefinitions are not expected to introduce new states, 
and hence we will only consider the case when r(— p^) has, in fact, no zeroes. Partially 
for simplicity, we also assume it to be analytic in the entire complex plane (except 
perhaps at oo). In this case one can write T{—p'^) as: 

r( V) = e^^"^'^ • (3.11) 
where 7(— p^) is any analytic function. 
3.2 Newtonian Potential and Asymptotic Freedom 

To derive the Newtonian potential one has to linearize the "generalized Einstein's field 
equations" corresponding to the action given in (L5) (see appendix |^ for details). In 



particular, the trace of the linearized field equations coupled to a point source reads: 

°° \ _ 1 

- 6 h = --□r(n)/i ~ 6{r) . (3.12) 




2 

Since we are only interested in cosmology, we only need to look at the potential for the 
scale factor or the spin piece, which is captured by the trace h of the linearized metric. 
Note that in the absence of any modification the massless propagator of the graviton can 
be captured alone by the □ in front of r(n). 

We now further make the assumption that the fields are varying slowly with time, so 
that □ — > V^, and after some massaging we obtain: 

Hr) ~ / 2^ ~ — / '^P ^'"^ ■ (3-13) 

J p^T{-p^) r p 

Note that the 1/r part arises mainly from the massless mode of the graviton. Moreover this 
is an important result of our paper as it tells us how higher curvature gravity in general 
modifies the propagator, and therefore the nature of gravity, as one spans very large to 
very small scales. For example it is clear that: 

• Any modification to the "Newtonian propagator" T that falls faster than 1/p can give 
rise to an asymptotically free theory of gravity. In particular, as long asj{—p'^) grows 
with large p^ , we will not only have a ghost free theory, hut also asymptotic freedom. 

For simplicity, let us restrict our attention to the case when 7 is a finite power series 
in w = —p^, and that there is a single scale m where the higher derivative effects become 
important: 

N /_ \ i 

^{w) = Y,h[^) , h^O{l) (3.14) 
i=i J 



Asymptotic freedom is guaranteed provided fcjv > 0. From Eqs. ( 3.13 ) and ( |3.14 ), it also 
follows that we recover Newtonian gravity in the large r limit, provided m > 10~^ eV, 
since we have tested gravity only at distances above the mm. scale, which corresponds to 
m ~ 10~^ eV. 
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3.3 Examples of Ghost and Asymptotically Free Theories 

For the purpose of illustration, let us choose the simplest example: 

'y{w) = -w^ r{w) = e-"" . (3.15) 

By construction this does not have any poles and is therefore ghost free^. 

Let us first find the action that corresponds to ( |3.15| ). From ( |3.6D , we can immediately 
match the coefficients by equating: 



which corresponds to 




6(i + l)! 

For the above choice the form of F{R) can be computed as: 

y ) |-1«D_D ^ D I ^ ^ 

i=0 

Next, let's look at the Newtonian limit. The Newtonian potential is given by 



(3.16) 



(3.17) 



F{R) = R - -RY } ' n'R = R--R[ \r. (3.18) 



h(r) = - / dp- e^P" = -erf ( - ) . (3.19) 

rj-^oo P r \roJ 

where rg is the ultraviolet cut-off, determined by the energy scale, governing the graviton 
exchange. 

Since lim j.^o6rf(r) ~ r and lim r^ooerf(r) = 1, we naturally obtain the kind of 
Newtonian potential we are seeking at large distances, i.e., 

lim h{r) -, (3.20) 

r— >oo T 

while for small r we are looking for an asymptotically free theory, with no singularity: In 
this case ^ 

lim h{r) ~ — = const. . (3-21) 

r^O ro 

This illustrates the weakness of gravity at short distances when r — > tq. 

Let us now illustrate how actions like ( |l.5D can also arise from symmetry consderations. 
For example, one realizes that a slight modification of the action corresponding to ( 3.15| ) 
leads to a form that is relatively simple in conformally fiat spaces. Consider introducing 
a scale factor into the action by a local (Weyl) scale transformation Qmn — ^ 't'^Qmrf 

This 

scalar (j) itself has weight 1 under a scale transformation. The combination (□ + ^R)4> also 
transforms covariantly, with weight 3. Then we can write 

S = jd'xV^ c/>3e'^-'(°+l«)0~2p + 1^)^ _ (3 22) 



®Such exponential propagators also appear in field theory duals of "discretised" string theory, for instance 
see [E9|. This type of modification to gravity has also been considered in pol]. 
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In the local scale gauge (p = 1, this becomes: 

Jd'xV^ e(°+i^) ^i?, (3.23) 

while for conformally flat spaces we can set the metric to be flat {(p itself being the scale 
factor) so 

f (fx (tf^e'^~'°(t>-'^U(t> , (3.24) 



where □ is now the flat one. It is easy to check that up to quadratic (in R) terms, ( |3.23| ) 
corresponds exactly to ( 3.18 ). 

4. Bouncing Cosmology from Asymptotically Free Gravity 



In the previous section we saw how nonperturbative gravity theories can encompass 
asymptotic freedom, which provides hope for resolving various singularities we encounter 
in the general theory of relativity. Here our main focus will be on the initial cosmological 
singularity. Although at first glance it may seem very difficult to understand the dynamics 
involving higher derivative actions such as ( |1.5| ), we will in fact be able to explicitly show 
approximate and exact bouncing solutions for a wide class of higher derivative gravity 
theory. 

First, we discuss why intuitively we expect to obtain bouncing solutions in asymp- 
totically free theories. Second, we provide exact bouncing solutions for ghost and asymp- 
totically free non-perturbative gravity actions described in ( |l.5D in the presence of a cos- 
mological constant. Although phenomenologically less appealing, it gives insight into the 
nature of gravity near the bounce. In the next section, we show how to obtain approximate 
bounce solutions in the absence of a cosmological constant, as well. We also highlight how 
to match the bouncing solutions to the standard radiation and matter dominated epoch. 

4.1 Intuitive Picture 

Intuitively, the reason we expect to have a bounce in theories with higher curvature gravity 
is because, although they do not have ghosts, gravity itself becomes weak at small scales. 
In the usual FRW cosmology, it is the attractive force of gravity which makes the universe 
contract, eventually to a singularity. However, now in the absence of strong gravity at 
short distances (or high energy densities) spacetime contracts but the internal pressure of 
matter can resist, therefore causing the universe to bounce from the contracting phase into 
an expanding one. 

In the Newtonian cosmology, the Hubble equation can be derived from the energy 
conservation equation of particle mechanics, i.e., K + V = E, 

lmr'^ + GN^r^pVir) = E. (4.1) 

Substituting V{r) = —1/r, we recover the analog of the Hubble equation: 



where the second term on the right corresponds to spatial curvature. 

Speciahzing to the case of a flat space (E = 0) and considering a general Newtonian 
potential arising from the higher curvature gravity action derived in previous section, we 
find 

H^^- = -^prV{r). (4.3) 

In an asymptotically free theory, at short distances the potential tends to a constant, so 
the right hand side of ( [4.3| ) goes to zero, suggesting — > 0, a condition signaling a bounce. 

One may wonder that perhaps one hits the singularity before the bounce can take 
place, because the gravitational energy vanishes only at r = 0. However, note that we have 
not included effects of pressure in ( [4.3[ ), which ensures that the bounce occurs at a finite 
scale factor. Unfortunately, Newtonian cosmology cannot capture the essence of pressure, 
but in the following subsections we will explicitly realize this intuitive picture in a large 
class of higher derivative actions. 

4.2 Hyperbolic Cosine Bounces 

To understand how one can deal with non-perturbative actions and equations of motion, 
let us look at the trace of the generalized Einstein equation of motion, given in appendix 
g, see O: 

oo / m \ 

G = G + GcoOR + ^ c„ en^+^i? + 2 X] n™-"i?n"i? = -4A , (4.4) 

m=l \ n=l / 

where A refers to the cosmological constant. Note that the trace equation involves D's 
acting on R. Thus the trick is to find a scale factor a{t) such that R{t) is a linear combi- 
nation of a set of functions Hi(t) with the property that the □ maps the functions Hi{t) 
onto themselves: 

OHi = MijHj , (4.5) 
where Mij is a constant matrix and the index i runs to some finite value. We can now find 



exact solutions, because all the terms in ( [4. 41) are at most quadratic in Hi and one has to 
only ensure that the overall coefficients of HiHj vanish^. 
A simple trial bounce solution is be given by^: 



a(t) = aocosh( J| t ] . (4.6) 



One finds the following: 



F = ,/| tanh(J|t) , (4.7) 



^This is not strictly necessary as sometimes they can be equated to matter sources with similar time 
dependences. 

*There exist other ansatze, such as a{t) oc e' , that can also provide non-singular bounce. For the 



purpose of illustration we shall stick to (4.6) 
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R 
\JR 



3w 



sech 



-{R + 3HR) 



w 

— t 



sech 



w ^ 



w[R — Gw] 



(4.8) 
(4.9) 



Naively one would have expected OR to contain terms proportional to sech^, but the contri- 
bution to such a term coming from the two pieces (the double derivative and the connection 
terms) cancel leaving only the sech^ piece, which is exactly what we wanted. Proceeding 
in a similar manner one finds 



In this example: 



w^'-^DR 



Ho = l, 



n / 0. 



sech 



(4.10) 



(4.11) 



Before we delve into bounces in non-perturbative gravity, let us briefiy mention what 
happens in the perturbative cases. Clearly, the prescription described above can be applied 
to the finite higher curvature theories as well to obtain hyperbolic cosine bounces. Albeit, 
one can show (see appendix ^ that they invariably involve negative energy density, in 
congruence with our earlier arguments. They either have ghosts in the gravity sector, or 
require the presence of ghost-like radiation with negative energy density. 

This also implies that we check that the radiation energy density corresponding to the 
bounce solution remains positive to ensure that the bounce is caused really by the higher 
derivative effects. This means that we need to study the Gqq equation, which is really the 
generalization of the Hubble equation in Einstein gravity. For an ansatz like (4.6) where 
the metric only depends on time, one finds (from (p.l|)): 



00 



p 1 °° 3 °° 

i^o^oo + --Fo.oo-aFo--J2 FnU^'R - - ^ FniU--^) . (4.12) 



n=l 



n=l 



After some some algebra we obtain 



00 



Lq + L2sech [ i — t ] + L^sech. \ — t\ + Lgsech 



,6 



w ^ 



(4.13) 



where one can conveniently express all the coefficients in terms of the kinetic operator 

r(n): 





w 
" 2"' 


L2{w) 


9w 


2 


La{w) 


9vu 






2 


Le{w) 


9w 




2 



1 — T{w) — icQW 



w 



3w 



2co + -T'{w) 





T{w) - 1 - T'{w)w 



w 



(4.14) 
(4.15) 

(4.16) 

(4.17) 
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Now, we wish to solve the 00 component of the generalized Einstein equation, ( C.l ), which 
reads: 



00 



Lq + L2sech I \ — i + L4sech' 



Too-^ 



- t ) + Lgsech" 



w 

2' 



A + posech { t 



(4.18) 



where in the 00 component of the energy momentum tensor we are compelled to include 
the cosmological constant term to match Lq, besides adding the radiation component. The 
radiation component must fall as sech^ following p ~ a~^. The above ( 4. IS] ) implies that 
there exists a solution if L2{w) = Lq{w) = or 1 — T{w) + T'{w)w = 0, and 1 — T{w) — 
4cqw — I = 0, or equivalently, 

T'(w) 



T{w)-T{0) 



w 



(r>.)> = ^r'(O) - i 



(4.19) 



(where ( ) is the average between and w). For the above solution, equating the sech^ 
term and the cosmological constant, we further find: 

PQ = ^(2wr'{w)-l) and ^ = ^- (4-20) 

To be a ghost free solution, we require: po > 

2wr'{w) - 1 > 0. (4.21) 

Thus, any kinetic operator F satisfying Eqs. ( [4.19| ) and ( [4.21 ) provides us with a ghost- 
free cosine hyperbolic bounce solution. We note that such cosine hyperbolic bouncing 
metrics and their generalizations (some of which also satisfy ( [4.5] )) were studied in [31]. In 
particular it was shown that although such metrics violate the dominant energy condition, 
causality is preserved and, as expected in bouncing universes, they can address the horizon 
problem. 

4.3 An Example 

In this subsection we provide a specific example which serves as an existential proof for a 
ghost free bounce. We argue that the above constraints, ( 4.19| ) and ( 4.21 ), can be satisfied 
for a wide range of F(t(;). Let us consider the following: 

j(w) = kiw - k2W^ + kiW^ , (4.22) 

where ki, k2, k^ are constants. Let us suppose that there exists a solution to ( 4.1S| ) at 
w = wq and 

T{wq) = Tq (4.23) 
Then the L2 and Lg constraints, along with ( 4.23| ), fix all the constants: 





1 


ki = 


2wq 




1 


k2 = 






^^0 - 


ki = 


1 


4 




^0 L 



(3Fo 



9Fg 



■2). 
4Fo 



Hi 

4F 



4Fo 
- 2 




2 In F 



InFn 



(4.24) 



- 14 - 



We compute the radiation energy density and we find: 



po = ^(2ro-3)>0^ro>^, (4.25) 



indicating that T{w) of the form ( 4.22| ), with coefficients given by ( 4.24| ), indeed gives rise 



to a ghost free bouncing solution, provided Tq > 3/2. wq remains a free parameter. 

The cosmological constant corresponding to this solution is given by A = We 
notice that we seem to require a large cosmological constant in order for the exact bounce 
solution to exist, just as in the example shown in appendix However, in the next 
section we will argue that a large non-vanishing cosmological constant is rather an artifact 
of our ansatz, i.e., a{t) oc cosh t, but is not a requirement for the bounce itself. 

5. Transition to FRW Universe 

It is clear that for large t the scale factor has a de Sitter solution: 

a(0-|e^/f^ (5.1) 

This corresponds to an inflation scenario that has no graceful exit. On the other hand, it 
seems a general feature of the exact "cosh" bounce solutions that they require a cosmo- 
logical constant (typically ~ 0{Mp)) for their existence. These two problems are related: 
Note as t ^ oo all the higher curvature terms vanish: 



n=0 ^ ' ^ 

and we are just left with Einstein's theory of gravity. The cosmological constant is pre- 
cisely required to sustain the inflation in this phase (J^^, and thereby ensures that the 
asymptotics of our ansatz is consistent. 

What is important is to realize that the cosmological constant is mainly required for 
the late-time consistency, see (^.1[), for our ansatz (4.6). It is certainly not necessary for 



having a bounce. In fact, we argue that its absence solves the graceful exit problem. 
5.1 Approximate Bounce Solutions 



We have seen how for large t all the higher curvature contributions in ( 4.12| ) vanish, leav- 



ing only the lowest order, Einstein-Hilbert term. In contrast, exactly at the bounce the 
contribution from the Einstein-Hilbert action is zero, while those coming from the higher 
order terms are finite. Thus, to understand the dynamics it seems natural to divide the 
evolution into two distinct regimes: 

1. near the bounce when the higher curvature terms dominate the evolution 

2. away from (both before and after) the bounce when the higher order terms can be 
ignored as compared to ordinary gravity and we have normal FRW evolution 
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In the next subsection we will study the transitions. Here we focus on how to obtain 
the approximate bounce solutions, i.e., phase 2. 

As compared to ( 4.18 ), the approximate equation that we now wish to solve contains 
only the higher derivative terms, 



w ^ 



(5.3) 



We neglected the contribution coming from ordinary Einstein gravity and assumed A = 0. 
In the next subsection we will check that it is indeed consistent to neglect the Einstein- 
Hilbert terms near the bounce. Now, it is clear from ( [4.14 -4.17) that Gqq^ is of the form: 



G§f^ = Lf ^sech^ [ \l — t] +^4sech^ ( a/tt * ) +^6sech° [ \l — t 



(5.4) 



where L4 and Lg are given as before by (4.16, 4.17), while 

9w'^ [1 - r(w) - icow 



tHD 

^2 



w 



(5.5) 



dU) then implies L^^{uj) = Lq{w) = 0, or l-T{w)+T' {w)w = and 1-T{w)-Acqw = 0, 
or equivalently, 

V(in\ - V(Vt\ 9 

(5.6) 



r»^£M^HrV))^^r'(o). 

w 6 



Moreover, equating with the stress energy tensor, we find: 

Po = ^:;-r [w) . 



(5.7) 




The constraints (|5.6| ) have simple geometrical interpretations. We can understand the 
solutions just by looking at the curve T{w) (see figure |l|). The function V{w) starts at 



r(0) = 1, and possibly passes through a finite number of minima and maxima on its way 
toward w ^ oo. The solution to the Lq constraint is given by the points at which tangents 
from the point {w = 0,r(0) = 1) touch the curve r(i(;). The constraint then simply 
says that the slope of this tangent has to be equal to two-thirds the slope of the curve at 

(y; = o,r(o) = i). 

Finally, the radiation energy density is also proportional to the slope of this tangent, 
so demanding a positive energy density solution simply means having a positive tangential 
slope. We note in passing that these constraints are very similar to the ones that we 
obtained for the exact bounce case. 

Again, as an existential proof one can check that for 

3 5 

ki = - ; k2 = --2ln2 and A:4 = l-ln2, (5.8) 



in (4.22), there exists a solution to the constraints at w = 1 and r(l) = 2, with positive 
radiation energy density po = 3/2. 

The pictorial understanding of the "bounce constraints" makes it easy to see what 
kind of T(w) can give rise to a bounce. However, we are still unable to clearly identify a 
connection of these constraints to the "necessity" or "sufficiency" of asymptotic freedom. 
For example, it is clear that T{w) corresponding to just the simple exponential as discussed 
in section |3.3| does not satisfy the bounce equations. The most likely reason is that we 
only looked at a very specific bounce solution, namely the hyperbolic cosine bounce, and 
therefore surely are missing more complicated bounces which may be present in some of 
the asymptotically free theories. We reserve a more detailed study of these issues to future 
research. 

5.2 Exiting the Bounce 



In this subsection we proceed as follows: First, we show why one can trust the approx- 
imate bounce solution that we have obtained in the previous subsection. We also estimate 
when the bounce solution breaks down and one expects a normal FRW universe to emerge. 
Second, we check that indeed away from the bounce the usual FRW cosmology is a good 
description and higher derivatives can be ignored. Finally, we obtain the full evolution 
of the universe by matching the scale factor and its derivatives at the boundary of the 
different regimes in a straightforward manner. For simplicity in the following we consider 
the action to be basically governed by a single scale as in ( p. 14 ) and therefore, unless there 
is some fine-tuning, one expects w ~ m?. Let us now investigate the different regimes of 
the evolution. 

As we argued earlier, (|4.6| ) is not an exact solution but it is still a good approximate 
solution at sufficiently high energy densities. One way to see this is to consider the left- 
hand side of (|5.3|) close to the bounce. In this regime, the higher order terms dominate 



over the Einstein-Hilbert term in the evolution equation (5^): Gqo = ^(1 — ^^^h^v^ *) 
while Gqq^ = ^^r'sech^y^ t. For example, exactly at the bounce {t = 0) Gqq vanishes 
while G^^ff w'^T' w is finite. 
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However, as t becomes large, all the higher order corrections get exponentially sup- 
pressed: 



sech^ (^^1 tj ~ e-Vf * , (5.9) 

while the zeroth order term coming from the ordinary Einstein-Hilbert action is not: There 
is a constant piece 



3w 



1 — sech^ I A I — t 



2 

Thus, the bounce solution breaks down approximately when 



|t|^oo 3w 



(5.10) 



Goo = G^o- ^ to = ±cosh-M^==_-^ , (5.11) 



and the Einstein-Hilbert term starts governing the evolution from here on. Since 
0(1), we have: 

- to ~ Oil) or to ~ \ - ~ (5.12) 
2 \ w 

Let us now check that indeed there exists late-time attractor solutions where higher 
derivative terms can be ignored and one has FRW cosmology. The transition from the 
bounce should lead to a radiation dominated epoch before and after the bounce. During 
these phases, we claim we have ordinary gravity coupled to an ideal gas of matter /radiation 
fluid satisfying 

p = ujp and p + 'iHp = ^ /) = /Jo , (5.13) 

Vflo / 

where — and + refer to the pre and post bounce solutions respectively. Therefore, at large 
times the solution reads: 

/ f \ 2/3(1+0;) / f \P 

a(t)^4{-) .„*(_) . (5.14) 

To see that this is a late-time attractor solution, let us compare Goo with G^ff. For the 
power law solutions, we obtain: 

Goo ~ i , (5.15) 

Thus, for large times |t| >> m~^, the higher order terms are negligible, and indeed ( 5.14| ) 
is a late-time attractor solution. However, when we reach |t| < m~^, the Gqq term (and in 
general Gq^) starts to dominate over the the Einstein-Hilbert term. The solution ( 5.14| ) 
is no longer valid. At this point we expect a transition to the approximate bounce solution 
given by ( |4.6| ). By matching the scale factor and its derivatives at the transition, we find: 

cosh { to ) , (5.17) 
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with to given by 



2 to tanh {y^ to 



(5.18) 



In particular, for radiation p = 1/2 and wM to ~ 0.77, consistent with (5.12). 



To summarize, the full evolution is given by 

cosh 



ait) 



cosh [y^ t 
cosh 



f to)(4^r fort<-to 

for to > t > —to 
for t > to 



(5.19) 



Note the scale is governed by the higher curvature action, while to is given by ( S.llj ). 
This completes our discussion on exiting an approximate bounce solution to the normal 
FRW universe. 



6. Future Directions 

In this work we have shown how non-perturbative effective actions of gravity are not 
only natural from the string theory perspective but also are able to realize a ghost and 
asymptotically free theory of gravity. Further, from our analysis it is clear that one can 
apply different techniques (like the Newtonian approximation, special ansatze, etc.) to 
understand the dynamics of the full non-perturbative theory. In particular, we obtained 
bouncing solutions and Newtonian potentials, which enabled us to argue how at short 
distances/times these theories can address singularities like the Big Bang singularity, while 
approximating ordinary General Relativity and FRW cosmology for large distances and 
times. 

Although we could only understand the transition from bounce to FRW cosmology 
in an approximate way, nevertheless this supports our view that bouncing cosmology can 
exist in a flat geometry without any requirements such as presence of ghosts or tachyonic 
states. 

In the light of observational data from WMAP, another important question would be 
the evolution of cosmological perturbations through the bounce, as well as their spectrum. 
It would also be interesting to investigate whether such a bounce leaves any distinctive 
signature in the production of gravitational waves, etc. Further note the fact that we could 
perform an exact analysis of non-perturbative actions provides us hope of trying to address 
other singularities, such as one encounters in the black hole physics. 

In brief, there are several interesting and important questions that we can now ask in 
the context of non-perturbative gravity. 
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A. Condition for a Singularity Free Universe in F{R) Cosmology 



We start by looking at scalar-tensor theories ( |2.5D , which can give rise to a monotonically 
expanding singularity free evolution, namely condition I. Note that the curvature scalar is 
given in terms of the scale factor, 




= 6 - + - . (A.l) 



From ( A.l| ) it is clear that we will have a singularity in the past if a = at a finite time. 
This can be seen very easily: Say a = at t = 0, then a ~ t*^, 5 > (note that a is 
analytic), therefore, this implies d ~ t^-^ and d ~ ^, which in turn implies R ^ t ^ and 
hence it diverges. 

Thus in a singularity free universe a{t) is bounded from below, but for condition I to 
hold we require d > at all times. This implies a{t) —> constant as t — > — oo, so that 
d, d — > as t ^ — oo. (See figure |2| for a possbile evolution of a{t).) 

The evolution of the scale factor reads, in the presence of a scalar field, 

H^ = t^ = \{K + V)^\p^, (A.2) 

^ = hv-2K), (A.3) 
a 6 

where K and V are the kinetic and potential energies of the scalar field. Since we argued 
that as t — > — oo, d, d — > 0, this then implies that /C, F — > as t ^ — oo, unless a ^ 0, too. 

Consider the first case, when the "particle" starts out with approximately zero kinetic 
and potential energies; then 

^ = -6H<i)^<0. (A.4) 
at 

This leads to being negative at all finite times, but this contradicts (|A.2| ). We are left 
with the only other possibility that a — > as t — > — oo such that d/a — > Hq ~ constant. 
This means 

a{t) *^ aoe-^«* . (A.5) 

In other words we have shown that an exponential inflation is the only way to realize 
the condition I. On the other hand it was shown in that such spaces are geodesically 
incomplete, and this is also the reason why it seems unlikely that such an evolution will 
be able to resolve the Big Bang singularity. However, they may be relevant for inflation. 
Therefore, let us proceed and try to pinpoint the kind of F{R) actions which can produce 
such an evolution. 

The evolution equations for the scale factor (A.2) and (A.3) suggest 



K and V constant > . (A.6) 

Specifically, since (/> — > 0, this means that — > 0, too. Then from the evolution equation 
for (p: 

4> + ^H^ = -^^^'^^0. (A.7) 
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This is the key result which tells us that to obtain an exponentially expanding solution we 
need dV/dcj) — > 0, a profile which asymptotes to a flat potential. This may happen only if 
dV/d(j) {) SiS (p ^ ±00. This criterion will dictate what kind of F{R) should be chosen. 
Let us first consider F to be a finite series: 

N 

F{R) = R + Y^ c„ii"+2 ^ (A.8) 

n=0 



with N > 2. Now, it is clear from (2.3) that in order for the mapping to work we need a 
sector of F{R) where F'{R) spans the interval [0, cxd], so that cj) spans (—00,00)^. 



F(R) 



F(R) 



Ro 
A 



B 



Ghost-free regimes in F(R) theories: For A, R > F^and for B, R < Rq 

V((^) 



a(t) 










t 



Possible evolution in a 
monotonically expanding Universe 



Potential corresponding to R gravity 



Figure 2: Cosmology with F{R) Theories 

There are two possible ways of realizing this criterion, which we have qualitatively 
sketched in Fig. 2. For concreteness let us restrict our attention to Fig. 2A. (Similar 
arguments can be made for Fig. 2B.) The relevant part of F{R) that can realize this 
starts at = i?o 7^ 0, which corresponds to = or ^ —00, and ends at i? — > cxd 
corresponding to F' ^ 00 or ^ 00. 

Note that we want to see whether V'{(j)) as (p 
to look at the behavior of 1^(0) as (p ^ ±00. Now as 
F{R) cnR^. The values of R and F{R) tend to: 



lb CXD. Therefore we really need 
00, ^ 00, and therefore. 



R 



1 



Nc 



1 

JV-l 



, V3(N-1) 



and F{R) 



CN 



1 



NCN 



N 
JV-l 



V2N 
, v^{JV-l) 



(A.9) 



®In reality the (j> field will be bounded by some cut-oflF for the validity of the low energy effective theory 
to work, but for the moment let us assume that (j) can take asymptotically large values. 
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After some algebra we obtain (see also 



1 



m^-^^^)""e-"^ with a = ^^^^<V2. (A.10) 
N \NcnJ VsiN-l) ^ ' 

Except for N = 2 this is an exponentially decaying function, i.e., V{(l)) — > as </> ^ oo. 
Thus this cannot work, as we required V{(l)) to asymptote to a positive constant ( |A.6| ). 
Next, let us look at (p ^ —oo or R ^ Rq, where 

F{R) « F{R^) + ^^^(i? - RoY , (A.ll) 
where is the first non-zero derivative of F at Rq. In this limit 

FiR)^Fm^^{^f'e^^. (A.12) 

In this case the asymptotic behavior of (j) depends on whether F{Rq) is nonzero or not. 
Transforming R, F{R), we obtain 

e if F{Ro) = 

As (j) ^ — oo, the potential grows exponentially large {V = ±00) rather than asymptoting 
to a constant as we want. Therefore, a generic F{R) cannot give rise to an exponential 
expansion or inflation. 

One exception arises (see ( A.l(]| )) when F{R) ^ i?^ as i? — > 00, i.e., 

F{R) =R + cqR^ - 2A, . (A.14) 

Indeed for such an action inflating solutions were found earlier |24]. Explicitly, one finds 
in this case 

For Co > and l/4co + 2A > this potential has the desired form to give rise to inflation. 

Even if one relaxes the condition that F has to be a finite series, the above conclusions 
essentially remain unchanged, as we will now show. First note that the argument about 
the behavior of V{(f)) as (p ^ —00 does not rely on the finiteness of the series, but only on 
the analyticity of F{R). Thus as — > —00, V{(j)) exponentially rises as before and cannot 
satisfy ( A. 6 ). We are left to consider what happens as — > 00. From the expression of 



the potential (12. q) , it is clear that if we want it to asymptote to a constant then at least 
i? — > e\/3 or F{R) — > e ^ (or faster). Consider the first case: This implies 

R > ^ = F'{R) =^ F{R) < R^ . (A.16) 
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For the second case 



F{R) > = F'^{R) ^ F{R) < . 



(A.17) 



To have an exponentially inflating universe it is necessary to have an F{R) which asymp- 
totes to ±R^ as i? — > ±00. Finally, we note in passing that adding radiation/matter does 
not change this conclusion. 

B. Cosine Hyperbolic Bounces in i? + cqR^ gravity 

Bounce solutions exist even when we consider only finite higher derivative corrections. 
However, when at least for one i > 0, q 7^ in ( |1.4| ), these theories necessarily have 
a ghost and therefore we do not consider them any further. Thus we are left with the 
combination F = R + cqR?. In this case the trace equation is given by (4.4): 



G 



Gnn + 2>Gy 



Qw — 2>w {1 — Qwcq) sech^ 



w 

— t 



-4A. 



The bounce solution is then given by: 

1 

w 



6co , and 4A = 



Co 



(B.l) 



(B.2) 



We however expect from the previous arguments that such a bounce solution is due to the 
"implicit" presence of ghost like radiation. To see this, one computes Gqc- In particular, 
one finds 



L4 = — 27cowsech 



9 



-sech 



t <0, 



(B.3) 



2 J 2 VV 2 

signaling a negative radiation energy density, or the presence of ghosts. 

We also note in passing that in the R^ type gravity, to get a bouncing solution, we 
require the presence of a typically large (corresponding to the scale of higher curvature 
corrections) cosmological constant. This is actually a generic feature of exact hyperbolic 



cosine bounces, and we addressed this issue in subsection 5.2. 



C. Newtonian Limit 



To obtain the weak field limit, we look at the higher curvature field equations directly. 

I F{R, BR, . . .) type of I 

Fg^u — Fo; tiu + QtiuOFo 



In [33 1 the field equations were derived for a general F{R, OR, . . .) type of Lagrangian: 

1 
2 



where 



Gnu = FnR 



^lu 1 



m=l 



dU^R 



(C.l) 



(C.2) 
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and as usual T^i, is the matter stress-energy tensor. It is easy to check that for ordinary 
Einstein gravity G^jy reduces to the usual Einstein tensor. Now, in order to derive the 
"Newtonian potential" corresponding to the scale factor it is sufficient to just look at the 
trace equation, which reads: 

oo 

G = g^'^Gf^u = FoR -2F + SDFo + 2 ^ FmW^R = -mS{r) . (C.3) 

m=l 

where we have coupled gravity to a point particle source of mass m. For F of the form 
( |1.5| ), one can show that: 

G = G + Y.^rn, (C.4) 

m 

where G = g'^'^Gfj,^ = —R, G^y being the usual Einstein tensor, 

Gq = QcqUR, (C.5) 
= Cm {qU^+^R + 2 a^-'^RWR^ V m / . (C.6) 

It is easy to check that g^y = rj^y is indeed the right vacuum solution of the full higher 
curvature equations (|C.3| ). Thus to obtain the Newtonian weak field limit we have to 
expand around the fiat space Minkowski metric, i.e., g^u = r]^y + h^jj] then 

G = -]pi, (C.7) 

where ^ 

h = hf'^; h^u = h^y - -'q^uh , (C.8) 

and we have used the Lorentz gauge h^y'^ = 0. 

It is a straightforward task to derive the potential for the scale factor,: The expression 
for G simplifies to 

~ \ - °° - 1 

G = — Uh + 3 ^ aU'+'^h = — UT{U)h = -m6{r) . (C.9) 
^ ^ 

Assuming h to be slowly varying with time, we find it to be given by: 

f „ p^P-^ 
h{r) ~ / (Tp 



p2r(— 

or 

- / dp— ^e'P"^ = — 

This completes our discussion on the Newtonian limit 
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